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^ i ABSTRACT 
Ph. 

^ ' Context. The mass assembly of galaxies can proceed through different physical processes. Here we report on the spectroscopic 

5_( ' identification of close physical pairs of galaxies at redshifts 2 < z < 4 and discuss the impact of major mergers in shaping galaxies at 

, these early cosmological times. 

^ . Aims. We aim to identify and characterize close physical pairs of galaxies destined to merge and use their properties to infer the 

contribution of merging processes to the early mass assembly of galaxies. 

Methods. We search for galaxy pairs with a transverse separation r,, < 25 h^'kpc and a velocity difference Av < 500 km s"' using 
the VIMOS VLT Deep Survey (VVDS) and early data from the VIMOS Ultra Deep Survey (VUDS) that comprise a sample of 1 1 1 1 
^ . galaxies with spectroscopic redshifts measurements at redshifts 1.8 < z < 4 in the COSMOS, ECDFS, and VVDS-02h fields. We 

' analyse their spectra and associated visible and near-infrared photometry to assess the main properties of merging galaxies that have 

' an average stellar mass Mt, = 2.3 x 10'" Mq at these redshifts. 

, Results. Using the 12 physical pairs found in our sample we obtain a first estimate of the merger fraction at these rdshifts, ~ 

. (15 - 20)%. The pair separations indicate that these pairs will merge within 1 Gyr, on average, with each producing a more massive 

galaxy by the time the cosmic star formation peaks at z ~ 1-2. From the average mass ratio between galaxies in the pairs, the stellar 
mass of the resulting galaxy after merging will be ~ 60% larger than the most massive galaxy in the pair before merging. We conclude 
that major merging of galaxy pairs is on-going at 2 < z < 4 and significantly contributing to the major assembly phase of galaxies at 
this early epoch. 

Key words. Galaxies: evolution - Galaxies: formation - Galaxies: high redshift - Galaxies: mergers 



- 1 — I 

X 

^ 1. Introduction 

" ■ ■ * Based on data obtained with the European Southern Observatory ^, ., . ^ , • , , 

Very Large Telescope, Paranal, Chile, under Large Programs 070.A- The contribution of different physical processes to galaxy mass 

9007, 177.A-0837 and 185.A-0791. Based on observations ob- assembly along cosmic time is still unknown, and a clear picture 

tained with MegaPrimeAlegaCam, a joint project of CFHT and describing how galaxies assemble, supported by observational 

CEA/DAPNIA, at the Canada-France-Hawaii Telescope (CFHT) evidence, has yet to emerge. Looking back at the average history 

which is operated by the National Research Council (NRC) of Canada, of a galaxy observed today, we are yet unable to identify how its 

the Institut National des Sciences de I'Univers of the Centre National stellai" mass has been acquired and which physical processes ai'e 

de la Recherche Scientifique (CNRS) of France, and the University possible contributors 

of Hawaii. This work is based in part on data products produced r,^, ri-- j 

, . J ^ .■ .5 . ,^ . . r The mass build-up of galaxies is expected to proceed 

at TERAPIX and the Canadian Astronomy Data Centre as part ot , , , . , „ , 

the Canada-France-Hawaii Telescope Legacy Survey, a collaborative trough a relatively small number of processes (for a summai-y, 

project of NRC and CNRS. This work is based on observations taken see e.g. Springel et al. 2005a; Silk & Mamon, 2012). New stars 

by the CANDELS Multi-Cycle Treasury Program with the NASA/ESA can form from the gas resei-voir of a galaxy, either acquired 

HST, which is operated by the Association of Universities for Research at birth or replenished from a more continuous accretion pro- 

in Astronomy, Inc., under NASA contract NAS5-26555. cess along the galaxy lifespan since formation. Major and minor 
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merging between galaxies is identified in numerous spectacu- 
lar examples in the local universe (e.g. in the RNGC catalogue, 
Sulentic & Tifft, 1973; Barnes and Hernquist, 1992). Merging 
is efficient at assembling mass, as it produces a significant in- 
crease in mass of up to a factor of two for equal mass mergers, 
for each merging event. Other processes are expected to mod- 
ulate the total mass gains. AGN and SNe feedback have been 
proposed as mechanisms capable of quenching star formation, 
as well as supporting winds capable at driving some mass frac- 
tion into the inter-galactic medium (IGM), hence reducing the 
increase in stellar mass from in-situ star formation (Silk, 1997; 
Murray et al. 2005; Cattaneo et al. 2009). The environment of 
galaxies is also expected to impact mass growth, with interac- 
tions between galaxies and the dense intra-cluster medium, like 
harassment or stripping, rapidly removing a significant part of 
the gas content of a galaxy (e.g. Moore et al. 1996). These pro- 
cesses are expected to ultimately combine along cosmic time 
to produce the mass distribution observed in the well-defined 
Hubble sequence of galaxy types in the nearby universe. 

In recent years, cold gas accretion fuelling star formation 
has received a focused attention, following numerical simula- 
tions (Keres et al. 2005). In this picture, cold gas flows along 
the filaments of the cosmic web into the main body of a galaxy 
to support vigorous star formation. This mechanism has been 
proposed as the main mode of galaxy assembly (Dekel et al. 
2009), and is often cited in the recent literature as the preferred 
scenario for galaxy assembly (e.g. Keres et al. 2009; Dijstra & 
Loeb, 2009; Bouche et al. 2010; di Matteo et al. 2012). However, 
as of today, only limited and indirect observational evidence ex- 
ist in support of this picture (Cresci et al. 2010; Kacprzak et al. 
2012), while detailed observational investigations have failed to 
identify direct supporting evidence (Steidel et al. 2010), which 
demonstrates the difficulty in directly identifying the accretion 
process at work. 

The merging of galaxies is another key process that con- 
tributes to galaxy assembly. The hierarchical growth of dark mat- 
ter halos is a key prediction of the KCDM model for galaxy 
formation (Davis et al. 1985; Springel et al. 2005b; Hopkins et 
al. 2006). In this picture, the merging of dark matter (DM) ha- 
los would not only lead to an increase of the DM halo masses, 
but also naturally lead to the merging of the galaxies associ- 
ated with each of the halos (Kauflmann, White & Guiderdoni 
1993). While galaxy assembly seems to produce the more mas- 
sive galaxies early in a seemingly anti-hierarchical downsizing 
pattern (De Lucia et al. 2006), it is nonetheless expected that 
merging of galaxies would continue to occur as DM halos con- 
tinue to merge along cosmic time. As the dynamical time-scale 
for halos to merge is of the order of 0.5-1 Gyr (Kitzbichler & 
White 2008; Lotz et al. 2010), it is expected that a massive halo 
today, as identified in DM halo merger trees from numerical N- 
body simulations, will have experienced several mergers since 
its formation. Documenting merging activity at difl'erent epochs 
can therefore shed light on the contribution of this process in 
assembling mass in galaxies. 

Evidence for merging is direct and well documented. 
Mergers have been identified since early days of photographic 
galaxy atlases and classified alongside the Hubble sequence of 
morphological types (e.g. the RC3 catalogue, de Vaucouleurs et 
al. 1991). The merging process of two disc galaxies has been 
proposed as one of the mechanism for producing early-type 
bulge-dominated galaxies, supported from simulations (e.g. 
Mihos & Hernquist 1996; Boumaud et al. 2005). Although this 
picture may be too simplistic (Bournaud et al. 201 1), it is, how- 
ever, clear that mergers do occur in the low redshift universe and 



that major mergers may lead to large modifications of the phys- 
ical properties of the galaxies involved. Mergers are identified 
either aposteriori morphologically from wisps, tails, or irregular 
shapes produced by on-going or post-merger dynamics, or apri- 
ori from the identification of pairs of physically bound galax- 
ies destined to coalesce. Going to high redshifts, using pairs in 
early merging stages, rather than post-merger remnants, remains 
the most robust way to derive a merger fraction (e.g. Le Fevre 
et al. 2000). This is because major merging pairs are easier to 
identify at these redshifts than the post-merger morphological 
signatures, which are of low surface brightness and may escape 
detection. The pair fraction can be transformed into a merger 
rate per volume or per galaxy (Patton et al. 2000; Kitzbichler 
& White 2008; Lin et al. 2008; de Ravel et al. 2009; Lopez- 
Sanjuan et al. 2011; Lopez-Sanjuan et al. 2013), using the dy- 
namical time-scale for a pair of galaxies with a given mass ratio 
and projected physical separation (Kitzbichler and White 2008; 
Lotz et al. 2010). The integrated merging rate over the lifespan 
of a galaxy since formation would then give the total amount 
of stellar mass assembled from the merger process. The secure 
identification of pairs and pair fractions at different redshifts is 
therefore an important observational measurement to perform. 
Spectroscopic redshift measurements of both members of the 
pair are required to eliminate the risk of background or fore- 
ground contamination along the line of sight and confirm that 
the pair is physically bound. 

At z ~ 1 pairs have been observed from deep galaxy spec- 
troscopic surveys (e.g. Lin et al. 2008; de Ravel et al. 2009). 
From the VIMOS VLT Deep Survey (VVDS), de Ravel et al. 
(2009) found that the merger fraction is higher by a factor ~ 3 
at z = 1 than in the local Universe and further showed that the 
merger rate and its evolution depends significantly on the stellar 
mass (luminosity) of the galaxy population. At redshifts beyond 
z ~ 1 only a few direct identifications of pairs and measurements 
of the merger fraction and merger rate exist. Lopez-Sanjuan et 
al. (2013) have reported a high pair fraction of ~ 20 - 22% at 
1 < z < 1.8 from 3D spectroscopy measurements in MASSIV 
(Mass Assembly Survey with Sinfoni in VVDS, Contini et al. 
2012). At higher redshifts Conselice et al. (2003) have used 
the CAS (concentration, asymmetry dumpiness) methodology, 
which relyies on image shapes and the expected signature of on- 
going or past mergers, to perform a measurement of the merger 
fraction up to z ~ 3. Cooke et al. (2010) have provided spectro- 
scopic identification of 5 pairs of galaxies in their LBG sample 
at z ~ 3, claiming that merging is triggering a significant part 
of the Lya emission. The number of confirmed pairs is therefore 
small beyond z ~ 2, and larger samples have to be assembled to 
enable measurements of the pair fraction and merger rate which 
is accurate to a few percent. 

Here we present a sample of galaxy pairs identified at z > 1.8 
in flie VUDS (VIMOS Ultra Deep Survey) and VVDS. The 
VUDS is an on-going survey (Le Fevre et al. 2013b, in prep.) 
with ultra-deep spectroscopy obtained with VIMOS on the VLT 
targeting galaxies with z > 2 in three well studied fields: 
file COSMOS, ECDFS and VVDS-02h (XMM-LSS/CFHTLS- 
Dl). The VVDS has been extensively discussed elsewhere (Le 
Fevre et al. 2005). We are using the "Final Data Release" of this 
survey as described in Le Fevre et al. (2013a, in prep.). These 
spectroscopic redshift survey samples are searched to identify a 
sample of real physical pairs with redshifts 1 .8 < z < 4, based on 
the observed projected separation Vp and velocity difference Av. 
We discuss the derived pair properties using all available spec- 
troscopy as well as visible and near-lR imaging and photometric 
data. The VUDS and VVDS spectroscopic redshift surveys are 
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described in Section|2] and the methodology to identify pairs and 
confirm that they are at close physical separation rather than a 
random projection along the line of sight is presented in Section 
[3] We then examine the pair properties in Section |4] We discuss 
the first detection of a significant number of major merging pairs 
at 1.8 < z < 4 and conclude in Section|5] 

Throughout this work, we adopt the "convergence cosmol- 
ogy" Hq = 100/z km 5-' Mpc-\ h = 0.7, Oo,a = 0.73 and 
f^o.Hi = 0.27. All magnitudes are given in the AB system. 

2. Spectroscopic observations 

To look for pairs, we have explored the VUDS, VVDS-Deep, 
and VVDS-UltraDeep surveys, providing a sample of galax- 
ies with spectroscopic redshifts measured with VIMOS on the 
ESO-VLT (Le Fevre et al. 2003). The VIMOS spectra have been 
obtained with 4.5h of integration for the VVDS-Deep survey, 
covering 5500 < A < 9350A, and 16h and 14h integrations in 
each of the LRBLUE and LRRED grism settings for the VVDS- 
UltraDeep and VUDS surveys respectively, covering a combined 
wavelength range 3600 < A < 9350A. The redshift accuracy 
with this setup is ^ 0.0005 x (1 -H z) (Le Fevre et al. 2013a). 

The VVDS-deep (Le Fevre et al. 2005; Le Fevre et al. 
2013a) and VVDS-UltraDeep (Le Fevre et al. 2013a) surveys 
are based on /-band magnitude selection with 17.5 < Iab ^ 24 
and 23 < Iab ^ 24.75, respectively. These two VVDS sur- 
veys are located in the VVDS-02h field centred at a2()00 - 
02/!26m00s and (52000 = -04deg30'00". 

The VUDS is an on-going spectroscopic survey also using 
VIMOS, targeting z > 2 galaxies in 3 fields: COSMOS, ECDFS 
and VVDS-02h. The baseline target selection for spectroscopy 
is using photometric-redshifts Zpiwi > 2.4 measured using all the 
photometry available in the survey fields using the Le Phare code 
(Ilbert et al. 2006). The photometric redshift accuracy obtained 
from the multi-wavelength data and calibrated on existing spec- 
troscopic redshifts is 0-37/(1+7) - 0.01 for magnitudes iab < 25 in 
the COSMOS field (see e.g. Ilbert et al. 2013). There is a slight 
degradation by a factor < 2 on the photometric redshift accuracy 
in the other two fields because of the lesser number of photo- 
metric bands observed. In addition, we supplemented the Zphoi 
selection by several color-selection criteria, adding those galax- 
ies likely to be in this redshift range, but not selected from the 
primary Zphot selection. Allowing for errors in Zphot, this selection 
provides a sample with 2 < z < 6, as described in Le Fevre et al. 
(2013b, in preparation). 

All VIMOS data have been processed with the VIPGI pack- 
age (Scodeggio et al. 2005). Following automated measure- 
ments, each galaxy has been examined visually and indepen- 
dently by two persons, each assigning a spectroscopic redshift. 
These measurements have been compared before assigning the 
final redshift measurement. A reliability flag has been assigned 
to each redshift measurement representing the probability for the 
redshift to be right. As consistently shown from the VVDS (Le 
Fevre et al. 2005), zCOSMOS (Lilly et al. 2007), and VIPERS 
(Guzzo et al. 2013) surveys, the reliability of flags reflect the 
statistical process of redshift assignment between independant 
observers and does not depend on the survey type or its intrin- 
sic quality, with flags 1, 2, 3, 4 and 9 having a probability to be 
right of ~ 50, 87, 98, 100, 90%, respectively (flags Ix with x=l, 
2, 3, 4, 9 indicate a broad line AGN; flag 2x are objects falling 
serendipitously in a slit next to a main target; and both have prob- 
ability distributions similar to the main flag categories). 

In addition to the VIMOS spectroscopic data, a large set 
of imaging data is available in the three fields covered by 



our pair search. The COSMOS field (Scovifle et al. 2007) has 
a full coverage with the HST/ACS F814W filter (Koekemoer 
et al. 2007) and includes, among other data, ugriz photom- 
etry from Subaru (Taniguchi et al. 2007), and more recently 
YJHK photometry from the Ultra Vista survey (McCracken et 
al. 2012). Spectroscopic redshifts from the zCOSMOS survey 
are also available (Lilly et al. 2007). The ECDFS is covered 
by the MUSYC survey in UBVRIz' (Gawiser et al. 2006) and 
partly by the CANDELS survey with the ACS and WFC3 on 
HST (Koekemoer et al. 2011). The VVDS-02h field has deep 
CFH12K BVRI photometry (Le Fevre et al. 2004), and even 
deeper CFHTLS ugriz photometry (e.g. Ilbert et al. 2006), as 
well as JHK photometry from the deep survey with CFHT- 
WIRCAM (WIRDS, Bielby et al. 2012). 



3. Pair identification 

Pairs have been identified using the projected transverse separa- 
tion rp and the velocity difference Av - c(vi - V2)/(l - V1V2) 
where v,, the normalized velocity of the x''' galaxy, is given as 
Vx = [(1 +Zv)^ - 1]/[(1 +Zv)^ + !]■ The survey samples have first 
been scanned for separation rp < 25h 'kpc and Av < 500 km 
s ' . We chose these separations as it is expected that pairs would 
merge in about 1 Gyr (e.g. Kitzbichler & White, 2008), meaning 
that a pair observed at z ~ 3 would have merged by the peak in 
star formation activity at z ~ 1.5-2 (e.g. Cucciati et al. 2012). 

Given the limitation of our ground based seeing observa- 
tions, we are not able to identify pairs superimposed along the 
line of sight or separated spatially by less than 1 arcsec, the av- 
erage image quality (FWHM) of the imaging data, correspond- 
ing to 5 h^'kpc at z ~ 3. In addition, pairs with larger rp or 
Av could also merge, albeit on a longer timescale, and with a 
lower probability (Kitzbichler & White, 2008). We will expand 
our pair sample to these larger separations and will correct for 
the geometrical limitations and selection function of the obser- 
vations, following e.g. de Ravel et al. (2009), when addressing 
the merger rate in a forthcoming paper 

We used all objects with reliability flags 2 to 9 (galaxies) and 
12 to 19 (broad-line AGN) for the primary targets, including a 
total of 1 1 1 1 galaxies, and an additional 811 galaxies with flags 
1 to 9 and 1 1 to 19, as well as 21 to 29 (objets falling serendipi- 
tously in the slit), for the companions. 

In looking for pairs, there is a possibility that the two objects 
identified could be two giant HII regions of the same galaxy, 
which would then be prominent in the UV rest-frame probed 
by the /-band, but which would appear as a single galaxy in the 
H or K band which probes redward of the D4000 or Balmer 
breaks at wavelengths less sensitive to contamination by younger 
stellar populations. On the contrary, the persistence of separate 
morphological components in the redder bands is an indication 
of the existence of two galaxies even if these are embedded in 
a diffuse light background. We therefore examined the /-band 
and ////T-band images of the pair candidates selected from their 
separation. As described in the next section, all of our pairs are 
well separated and both galaxies are seen from bluer to redder 
bands. We are therefore confident that we are dealing with true 
physical pairs rather than double HII regions of the same galaxy. 

In this process, we have identified 12 pairs with redshifts 
from z - 1.82 to z = 3.65. All primary and companion galaxies 
have flags> 2 except one galaxy (ID:520478087) with a lower 
reliability flag 1 . Our pair sample therefore benefits from reliable 
spectroscopic identification. The pair properties are discussed in 
the next section. 
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4. Pair properties 

From the multi-wavelength dataset, we have derived the main 
properties of the galaxies in the identified pairs. The stellar mass 
of each galaxy of a pair has been derived from SED fitting of 
the available photometric data at the measured redshift (see e.g. 
Ilbert et al. 2010) and range from 10'' to 10" Mq (TablelD, with 
uncertainies by up to a factor 2 (e.g. Bolzonella et al. 2009) de- 
pending on the number of bands used, or the IMF. The stellar 
mass ratio are used to separate "major-merging pairs", with a 
mass ratio between the two galaxies of 1 < Mi /M2 < 4, from 
"minor-merging pairs", for which Mi /M2 > 4. We have checked 
that the SED-derived mass ratio for each pair is consistent with 
the TiT-band or //-band flux ratio of the two galaxies, as these 
fluxes probing above Aresi k, 4000A can be considered as a rough 
proxy for stellar mass (Bruzual and Chariot, 2003), and we have 
used the error on the /T-band or //-band flux ratio as a proxy 
for mass ratio errors. We find that the mass ratio of the major 
merging pairs are in the range 1.2 - 2.3 and therefore, given 
the Icr uncertainties on stellar mass ratio estimates reported in 
Table |2] none of the major merger pairs could be misclassified 
minor mergers. Taking into account the mass ratio between the 
two galaxies in a pair, and the separations Av and r^, we use the 
formalism of Kitzbichler and White (2008, formula 10) to com- 
pute the merging time-scale for each pair T^erg, as described 
in de Ravel et al. (2009). These timescales are compatible with 
the Lotz et al. (2010) estimates, as discussed in Lopez-Sanjuan 
(2011). We list the main pair properties in Table |2] including 
Tmerg and the redshift Zassembfy by which the pair would have 
merged into a single galaxy, taking into account the redshift of 
the pair and T„,^rg- Images and spectra of each pair are presented 
in Figures[T]to[T2] We describe the properties of each pair below. 

COSMOS-511001467 A/B (Figure[D: Redshifts of the two 
galaxies in this pair have been measured from two different 
VUDS observations at zi = 2.0970 and zi = 2.0961, for a ve- 
locity difference of Av - 87 km s"'. The two galaxies are sepa- 
rated by 18.4 /j^'kpc, and are easily visible both in the HST/ACS 
F814W image and in the Ultra Vista near-infrared images. The 
most massive galaxy has a stellar mass of 0.3 x 1O'°M0 and 
the stellar mass ratio between the two galaxies is estimated to 
be M1/M2 - 2, hence a major merger The South- West galaxy 
shows a compact nucleus surrounded by a faint difiiise compo- 
nent, while the North-Eastem galaxy is seen as two components 
in both HST/ACS F814W and Ultra Vista images. The specti-a of 
the two galaxies have a UV flux rising to the blue up to Lye, typi- 
cal of star forming galaxies at these redshifts. Given the physical 
separation and mass difference, and following the prescription 
of Kitzbichler & White (2008), this pair is expected to merge 
within 2.3 Gyr. 

COSMOS-510788270 A/B (Figure|2): Two galaxies sepa- 
rated by 15.4 /1 'kpc are at a redshift z - 2.9629, with a velocity 
difference Av = 38 km s"', measured from two different VUDS 
observations. The most massive galaxy has a mass O.7x1O'"M0 
and the mass ratio between the two galaxies is estimated to be 
M1/M2 = 1.5, a major merger One of the spectra shows Lya 
in emission with EW(Lya)rest - 15 A, while Lyo- is observed in 
absorption for the other galaxy. Images are barely resolved at the 
HST/ACS resolution, showing a slightly extended and elliptical 
light distribution. This pair is expected to merge within 1.5 Gyr 

COSMOS-510175610/510778438 (Figure |3): This pair at 
a redshift zi = 3.0939 is separated by 9.4/!"'kpc and Av = 161 
km s"^ measured in the same slit from VUDS observations. The 
brightest and most massive galaxy has a compact, although ir- 
regular, morphology, while the companion to the West, which is 



confirmed at the same redshift, is more diffuse and has a low sur- 
face brightness. Another companion is visible about 3/1 'kpc to 
the North-East, and although its photometry is compatible with 
the redshift of the pair, no spectroscopic redshift information is 
available to confirm that it is physically linked to this system. 
Both spectra have been obtained from the VUDS survey. One of 
the VIMOS spectra shows a weak Lya in emission, while the 
other is purely in absorption. The mass of the brightest galaxy 
is 5 X IO'^Mq, and the mass ratio is M1/M2 = 6.3, considered 
to be a minor merger This pair is expected to merge within 0.7 
Gyr 

VVDS-02h-520452183/520450423(Figure|4ll: The spectra 
of these two galaxies obtained from the same VVDS slit show a 
broad-line AGN with z - 1.8370 and an absorption-line galaxy 
with z - 1 .8333, for a velocity difference Av = 391 km s"' . This 
pair is separated by 22.3 /i^'kpc, and the estimated mass ratio 
is Ml IM2 - 1 .9, with the most massive galaxy, the AGN host, 
having a stellar mass 5.6 x 10"*Mq. This estimate is, however, 
quite uncertain given the presence of the AGN. Another object is 
observed in between the two main galaxies but no redshift infor- 
mation is available. The AGN host is compact, barely resolved 
at the seeing of the best CFHTLS image {FWHM = 0.6 arcsec), 
while the companion galaxy is slightly elongated and irregular 
in shape. Given the physical separation and mass difference, this 
pair is expected to merge within 1.1 Gyr, although this estimate 
is affected by the uncertain mass estimate of the AGN host. 

VVDS-02h-520478238/520478087 (Figure |5): The two 

galaxies observed spectroscopically in the same VUDS slit have 
Zi = 2.2460, Z2 = 2.2463 for a velocity difference Av - 28 
km s"', and are separated by = 25 /i^'kpc. They are part of 
a group of four galaxies within 5 arcsec (30/i"'kpc), with pho- 
tometric redshifts compatible with the redshift of the pair. The 
eastern component is the most massive with = 2.0x10'°Mq. 
The mass ratio. Mi /M2 = 2, indicates a major merger and results 
in an expected time to merge within 1.8 Gyr. 

VVDS-02h-910260902/910261083 (Figure |6): The two 

galaxies, spectroscopically observed in the same VVDS slit with 
zi = 2.3594 and Z2 = 2.3588, are identified both in the CFHTLS 
/-band image and the /T-band WIRDS image, with rp - 12.3 
/j^'kpc and Av - 54 km s"' separation. The slit was placed on 
the centroid of the blended CFHTLS image of the two galaxies, 
but still included a significant fraction of the flux of both galax- 
ies to yield two well separated spectra. The South-East galaxy 
is the most massive with 7.9 x 10'"Mq and the mass ratio is 
M1/M2 - 2.4, the later indicating a major merger While the 
massive galaxy shows a symmetric elongated shape, the North- 
West galaxy shows an irregular morphology, and the two are 
connected by a faint bridge. Both galaxies show Lya in emission, 
with EW{Lya)rest - 50A and 60A, and integrated Lya line flux 
Liya = 10"^^ and 3 x 10^^ erg s"' respectively, indicating strong 
star formation at the level of 1 - 2MQ/yr (see e.g. Cassata et al. 
2011). Given the observed separation of this pair, it is expected 
that it will merge within 0.6 Gyr 

VVDS-02h-910302317/910302929 (Figure E): The two 

main galaxies observed in spectroscopy in the same VVDS slit 
are respectively at zi = 1.8171, and Z2 - 1.8154, separated by 
8.7 /z^'kpc along the East-West direction, and by Av - 181 km 
s"'. The brightest / most massive galaxy has amass 1.8x10"'Mq 
and is made of two main components, well visible both in the i- 
band and /T-band images. The mass ratio between the two galax- 
ies in the pair is M1/M2 - 1.7, which indicates a major merger 
This pair is expected to merge within 0.6 Gyr 
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ECDFS-530034527 (Figure|8): This is the highest redshift 
for which we have obtained a spectroscopic confirmation of a 
physical pair, with both galaxies at the same redshift z - 3.6500, 
separated by 6.8 /i"' kpc, as measured from the same VUDS slit. 
In Figure |8] another companion is obsei-ved to the North-East, 
but there is no confirmation of its redshift. The stellar mass of 
the brightest and most massive galaxy is 1.2 x 1O'"M0, and the 
mass ratio with the companion is Mi /M2 = 1.3, which indicates 
a major merger. The two spectra obtained in the VUDS survey 
show absorption line spectra, with only a weak Lya emission for 
the brighter galaxy. This pair is expected to merge within 0.5 
Gyr. 

ECDFS-530050663 (Figureg): This pair is made of 2 com- 
ponents identified in the //-band CANDELS image, with red- 
shifts zi = 2.3250 and zi = 2.3244 measured in the same VUDS 
slit. The stellar mass of the main component is 0.4 x 1O'"M0, 
and the mass ratio between components is M1/M2 = 2: a ma- 
jor merger. A third component is obseiTed to the Noith-East, 
all three components being embedded in a low surface bright- 
ness emission. The pair identification is therefore ambiguous, as 
this configuration might be indeed a merger at an advanced stage 
surounded by tidal debris, or the result of three giant star forming 
regions in a single forming galaxy. Under the merger hypothesis, 
the two main components are expected to merge within 0.6 Gyr. 

ECDFS-530042814/530042840 (Figure Hfi): The two 
galaxies are well-visible in the composite BVR MUSYC im- 
ages and the /T-band image. These galaxies are at a mean red- 
shift z - 2.9903 with a velocity difference Av - 278 km s"' as 
measured in the same VUDS slit. These two main galaxies are 
separated by 6 /; 'kpc, and a third component is observed 6.5 
/i^'kpc to the North-East. While the main galaxy shows Lya in 
absorption, the other galaxy measured with VIMOS shows Lya 
in emission, both spectra having a UV slope indicating strong 
star formation. The mass of the main galaxy is estimated to be 
0.12 X 1O'"M0, and the mass ratio is M1/M2 - 1.2 almost an 
equal mass major merger. This pair is expected to merge withins 
1.0 Gyr. 

ECDFS-530046916 (Figure HD: Two close but well sepa- 
rated images of two galaxies are observed in both HST z-band 
(F850W) and //-band (F160W) images, with a redshift around 
z - 2.868 measured in the same VUDS slit. The brightest / most 
massive galaxy has a regular morphology compatible with a disc, 
while the companion galaxy has an irregular shape. The mass of 
the main galaxy is estimated to be 1.4 x 10"^ Mq, and the mass 
ratio is M1/M2 = 2.3, indicating a major merger. This pair will 
probably merge within 0.5 Gyr. 

ECDFS-530036900/26931 (FigureHSll: Two well sepai-ated 
galaxies are observed in the HST z-band (F850W) and //-band 
(F160W) images from CANDELS, with 7.0 /i^'kpc and Av = 78 
km s"' separation at a redshift of z\ = 3.3300. One galaxy is 
quite compact, while the companion shows a sharp point-like 
component surrounded by a nebulous extension. The compact 
galaxy has been obsei-ved by VUDS and shows strong Lya emis- 
sion, while the fainter galaxy has been observed as part of the 
VIMOS-GOODS survey (Popesso et al. 2009). The mass of the 
main galaxy is estimated to be 1.3 x IO'^Mq, and the mass ratio 
is M1/M2 - 1.2, an almost equal mass major merger. This pair 
will probably merge in less than 0.6 Gyr. 

5. Discussion and conclusions: major merging 
pairs at 2 < z < 4 

We have identified 12 pairs of galaxies with redshifts 1.81 < 
z < 3.65 from the on-going VUDS and VVDS spectroscopic 



surveys. Both components of the pairs have a confirmed spec- 
troscopic redshift obtained with VIMOS on the VLT and, there- 
fore, comprise a unique sample of true physical pairs at these 
redshifts. The galaxies in our sample span a large mass range 
from 10^ 'Mq to 10"Mq. The majority of galaxies in our pairs 
show signs of strong star formation, a common property at these 
redshifts, as evidenced by either strong Lya emission, strong UV 
continuum, or both, and with one of them showing AGN activ- 
ity. The mass ratio of the merger is in the range 1 < Mi /M2 < 6, 
with 11 of 12 pairs satisfying a major merger pair criterion 
1 < M1/M2 < 4. The VUDS sui-vey is still going-on, and deriv- 
ing the survey selection function is beyond the scope of this pa- 
per, however, we use these pairs to derive a preliminary estimate 
of the major merger pair fraction. After correcting for the s elec- 
tion function and geometrical effects we derive a first tentative 
measurement for the major merger fraction of /„j,r^ ~ (15-20)% 
at 1.8 < z < 4.0. This value is comparable to the one observed 
from pair identification at 1 < z < 1.8 in the MASSIV survey 
(Lopez-Sanjuan et al. 2013), and measurements obtained from 
CAS morphological analysis up to z = 3 (Conselice et al. 2003). 
Our merger fraction estimate will be further refined as the VUDS 
survey proceeds. 

Using the observed projected spatial and velocity separa- 
tions, along with the prescription from Kitzbichler and White 
(2008), we have computed the merging timescale for each pair 
(Table |2]l. With separations from 6 to 25/i^'kpc and mass ra- 
tio 1/6 to 1, the average merging timescale of our sample is 
< Ti„eig >= 1 Gyr and the median is 0.7 Gyr. Most of these 
pairs will therefore have merged before the peak of star forma- 
tion at z ^ 1 .5 (e.g. Cucciati et al. 2012). Given the average mass 
ratio of 1 .75 for the 1 1 major merger pairs, the main galaxies in- 
volved in these mergers will have increased their stellar mass by 
~ 60% from z ~ 3 to z ~ 1.5 from the merging process alone. 

Our obsei-vations therefore provide unambiguous evidence 
of major merging occurring at 2 < z < 4. It is clear that hi- 
erarchical assembly, with massive galaxies being built from the 
merging of less massive ones, is at work at these redshifts and 
that major merging is contributing to the assembly of mass in 
galaxies at early times. This mass assembly simply results from 
the sum of the masses in each galaxy in a merging pair, a simple 
and effective way to increase mass at each merging event. An 
additional increase of stellar mass from star formation triggered 
by the merging process is also possible, with a range of mass 
production identified in the literature, ranging from relatively 
large (e.g. Kocevski et al. 201 1) to more limited star bursts (e.g. 
Mullaney et al. 2012), depending on the duration and strength 
of the merger induced burst. Merging is therefore a clear path 
to move low mass galaxies towards the higher end of the mass 
function, which contributes to the evolution of the stellar mass 
function (e.g. Ilbert et al. 2012; Ilbert et al. 2013). As identified 
by de Ravel et al. (2009) and Lopez-Sanjuan et al. (2013), galax- 
ies in the local universe with a mass 10'°'^Mq will have assem- 
bled 30 - 40% of their mass from major merging since z ~ 1.5. 
Our results demonstrate that merging was also contributing to 
mass assembly in galaxies at even earlier times. 

These results are to be placed in the context of the currently 
favored picture of galaxy assembly, with cold accretion playing a 
key role in building-up mass in galaxies (e.g. Dekel et al. 2009). 
While cold accretion is observed in simulations as the primary 
mechanism fueling galaxies with enough gas mass to sustain the 
strong, continuous star formation necessary to explain the star 
formation rate history, the observational evidence for cold ac- 
cretion at these redshifts remains elusive. While our results con- 
tribute to bring evidence that merging is at work, it is important 
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Table 1. Pair properties 



Field 


Pair Galaxies 
ID numbers 


Redshift Stellar Masses'" 


Mass ratio'^' 


Flux ratio'^' 


' l> 

h^'kpc 


Av<5' 
km s" 


7^(6) 
1 fnerg 

Gyr 


assembly 


COSMOS 


511001467a&b 
510788270a&b 
510175610/510778438 


2.0970 
2.9629 


0.3/0.15 
0.5/0.7 

OS/SO 


2.0 
1.4 
6.3 


1.9 ±0.2 
1.7 ±0.3 
4.2 ± 0.2 


18.4 ±0.3 
15.4 ±0.2 
9.4 ±0.2 


87 
38 
161 


2.3 
1.5 
0.7 


1.1 
1.8 
2.3 


VVDS-02h 


520452183/520450423 
520478238/520478087 
910260902/910261083 
910302317/910302929 


1.8370 
2.2460 
2.3594 
1.8171 


5.6/3.0 
2.0/1.0 
3.4/7.9 
1.1/1.8 


1.9 
2.0 
2.3 
1.6 


2.2 ±0.1 
1.8 ±0.4 
2.0 ±0.1 
1.4 ±0.1 


22.3 ± 1.2 
25.3 ± 1.2 
12.3 ± 1.1 
8.7 ± 1.2 


390 
28 
54 
181 


1.1 
1.8 
0.6 
0.6 


1.3 
1.3 
1.9 
1.5 


ECDFS 


530034527a&b 

530050663a&b 

530042814/530042840 

530046916ac&b 

530036900a&b 


3.6500 
2.3250 
2.9903 
2.8684 
3.3300 


1.2/0.9 
0.4/0.2 
0.12/0.1 
1.4/0.6 
1.3/1.1 


1.3 
2.0 
1.2 
2.3 
1.2 


1.5 ±0.1 
1.9 ±0.2 
1.3±0.1 
2.5 ±0.1 
1.5 ±0.1 


6.6 ± 0.2 
5.5 ±0.3 
6.0 ±0.2 
6.5 ±0.2 
7.0 ±0.2 



54 
278 
279 
76 


0.5 
0.6 
1.0 
0.5 
0.6 


2.9 
1.9 
2.1 
2.4 
2.6 



Table 2. Stellar mass estimates from SED fitting; Stellar mass ratio between the two galaxies in the pair, errors on the mass 
ratio have been computed from the error on the ratio of the H or K band flux measurements (see text); Flux ratio as measured on 
H or K-band r^: transverse separation between the two galaxies in the pair, errors in have been conservatively estimated taking 
a one pixel error on the difference in the centroid measurement of each galaxy in the pair; Av: velocity separation along the line 
of sight; errors in Av are estimated to be < 100 km s"' when the two redshifts are from the same slit (10 pairs), and ~ 200 km s"' 
when redshifts are from two different slits (2 pairs); T„erg- timescale for the pair to merge, using the Kitzbichler & White (2008) 
prescription; from errors in mass ratio, and Av errors on r,„f,-g are ~ 10 - 20%; Zassembh'- redshift by which the two galaxies 
will have merged, obtained combining the observed redshift and the Az corresponding to the merger timescale. 




Fig. 1. Pair COSMOS-5n001467 A/B: 5" x 5"sum of UltraVista YJHK images (top-left) and HST F814W image (bottom-left). 
The location of the 1 arcsec width VLT/VIMOS slits is shown by the rectangles. The more massive object is labeled as A. Right 
Panel: VIMOS spectra for both components in the pair The spectra have been arbitrarly shifted in flux to avoid overlap. 



to note that merging is not exclusive to other mass assembly pro- 
cesses. The high gas fraction of high redshift galaxies (Tacconi 
et al. 2010; Daddi et al. 2010) implies that other processes must 
also be at work, such as the continuous formation of stars po- 
tentially fueled by cold accretion. The total mass growth of a 



galaxy must therefore be coming from all these different pro- 
cesses at work. By precisely knowing the contribution of merg- 
ing to global mass assembly, measured, e.g., by the growth of the 
stellar mass density, one would ultimately be able to place upper 
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Fig. 2. Pair COSMOS-5 10788270 A/B: 5" x 5"sum of Ultra Vista YJHK images (top-left) and HST F8 14W image (bottom-left). 
The location of the 1 arcsec width VLT/VIMOS slits is shown by the rectangles. The more massive object is labeled as A. Right 
Panel: VIMOS spectra for both components in the pair. The spectra have been arbitrarly shifted in flux to avoid overlap. 
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Fig. 3. Pair COSMOS-5 10175610/5 10778438: 5" x5"sum of UltraVista YJHK images (top-left) and HST F814W image (bottom- 
left). The location of the 1 arcsec width VLT/VIMOS slit is shown by the rectangle. The more massive object is labeled as A. Right 
Panel: VIMOS spectra for both components in the pair. The spectra have been arbitrarly shifted in flux to avoid overlap. 



limits to the total mass growth from other processes, including 
cold gas accretion. 

The integrated contribution of merging processes to the com- 
plete history of mass assembly requires knowledge of the evolu- 
tion of the merger rate since early times (i.e., significantly be- 
yond z ~ 2). Building on the sample presented in this paper, the 
VUDS survey, when complete, will enable a robust measure- 



ment of the merging rate out to z ~ 4, and an estimate of the 
total amount of mass assembled by the merging process since 
the early universe. 
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Fig. 4. Pair VVDS-2h-520452183/520450423: 6.5" x 6.5"composite WIRDS JHKs image (top-left) and CFHTLS i-band image 
(bottom-left). The location of the 1 arcsec width VLT/VIMOS slit is shown by the rectangle. The more massive object is labeled as 
A. Right panel: VIMOS spectra for both components in the pair. The spectra have been arbitrarly shifted in flux to avoid overlap. 
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Fig. 5. Pair VVDS-2h-520478238/520478087: 6.5" x 6.5"composite WIRDS JHKs image (top-left) and CFHTLS i-band image 
(bottom-left). The location of the 1 arcsec width VLT/VIMOS slit is shown by the rectangle. The more massive object is labeled as 
A. Right panel: VIMOS spectra for both components in the pair. The spectra have been arbitrarly shifted in flux to avoid overlap. 
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A. Right panel: VIMOS spectra for both components in the pair. The spectra have been arbitrarly shifted in flux to avoid overlap. 
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Fig. 8. Pair ECDFS-530034527 A/B: 5" x 5"HST/WFC3 F160W CANCELS image (top-left) and composite B VR image from the 
MUSYC survey MUSYC (bottom-left). The location of the 1 arcsec width VLT/VIMOS slit is shown by the rectangle. The more 
massive object is labeled as A. Right panel: VIMOS spectra for both components in the pair. The spectra have been arbitrarly shifted 
in flux to avoid overlap. 
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survey (bottom-left). The location of the 1 arcsec width VLT/VIMOS slit is shown by the rectangle. The more massive object is 
labeled as A. Right panel: VIMOS spectra for both components in the pair. The spectra have been arbitrarly shifted in flux to avoid 
overlap. 
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Fig. 10. Pair ECDFS-5300428 14/530042840: 5" x 5"composite JK image (top-left) and composite BVR image from the MUSYC 
survey (bottom-left). The location of the 1 arcsec width VLT/VIMOS slit is shown by the rectangle. The more massive object is 
labeled as A. Right panel: VIMOS spectra for both components in the pair. The spectra have been arbitrarly shifted in flux to avoid 
overlap. 
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Fig. 11. Pair ECDFS-530046916 A/B: 5" x 5"HST/WFC3 F160W image (top-left) and HST/ACS F850W image from the 
CANCELS survey (bottom-left). The location of the 1 arcsec width VLT/VIMOS slit is shown by the rectangle. The more massive 
object is labeled as A. Right panel: VIMOS spectra for both components in the pair. The spectra have been arbitrarly shifted in flux 
to avoid overlap. 
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Fig. 12. Pair ECDFS-530036900 A/B: 5" x 5"HST/WFC3 F160W image (top-left) and HST/ACS F850W image from the 
CANCELS survey (bottom-left). The location of the 1 arcsec width VLT/VIMOS slit is shown by the rectangle. The more massive 
object is labeled as A. Right panel: VIMOS VUDS spectrum for component A. The redshift of component B has been measured by 
Popesso et al. (2009) and not shown in the figure. 
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